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C , Abstract. A new method is developed for estimating photometric redshifts, using re- 
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alistic template SEDs, extending over four decades in wavelength (i.e. from 0.05 fim to 

00 

1 mm). The template SEDs are constructed for four different spectral types of galax- 
I ies (elliptical, spiral, irregular and starburst), satisfying the following characteristics: a). 
, they are normalised to produce the observed colours of galaxies at z ~ 0; b). incorpo- 
rate the chemo-photometric spectral evolution of galaxies of different types, in agreement 
2 with the observations; c). allow treatment of dust contribution and its evolution with red- 
^ shift, consistent with the spectral evolution model; d). include absorption and re-emission 
^ [ of radiation by dust and hence, realistic estimates of the far-infrared radiation; e). in- 
. elude correction for inter-galactic absorption by Lyman continuum and Lyman forest. 
Using these template SEDs, the photometric redshifts are estimated to an accuracy of 
Az = 0.11. 

The simultaneous and self-consistent modelling of both the photometric and chemical 
evolution of galaxies and the effect of dust, makes this technique particularly useful for 
high redshift galaxies. The effect on the estimated photometric redshifts, due to assump- 
tions in the evolutionary population synthesis models, are investigated and discussed. 
Also, the degeneracy in the predicted photometric redshifts and spectral types are ex- 
amined, using a simulated galaxy catalogue. 
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1. Introduction 

Recent multi-waveband galaxy surveys, carried out from UV to radio wavelengths, have 
identified a population of rapidly star forming galaxies in the range < z <3 (Sullivan et 
al 1999, Lilly et al 1996, Cowie et al 1997, Rowan-Robinson et al 1997, Barger et al. 1998; 
Hughes et al 1998, Mobasher et al 1999). In particular, these studies confirm a relatively 
higher rate of star formation in the past, as supported by the discovery of a population 
of massive, starforming galaxies at 2; 3 (Madau et al 1996; Steidel et al 1996). These 
objects are likely to be progenitors of the present day galaxies (Giavalisco et al 1996; 
Lowenthal et al 1997) and hence, a statistical study of this population from the epoch 
of galaxy formation to the present, gives clues towards scenarios of the formation and 
evolution of galaxies (Fukugita et al 1996). This can also constrain the star formation 
history of galaxies out to 2: ~ 3. Such studies require redshift information for a population 
of starforming galaxies at faint levels. Recently, the depth of the available surveys is 
greatly extended by the Hubble Space Telescope (HST) observations of the Hubble Deep 
Fields (HDF). The large spectral coverage of the HDFs provide a unique opportunity to 
study evolutionary properties of faint galaxies. 

Ground-based spectroscopic measurements of the brighter (/ < 25 mag.) sub-sample 
of the HDF have been performed (Cohen et al 1996; Steidel et al 1996; Lowenthal et 
al 1997; Zepf et al 1996). However, for the fainter galaxies in the HDF, redshift mea- 
surements are more difficult with spectroscopic features almost impossible to identify. 
For these objects, the photometric redshift technique (Loh and Spillar 1986; Connolly 
et al 1995) is faster than its spectroscopic counterpart and applicable to much fainter 
magnitudes. This is due to a larger bin size in photometry compared to spectroscopy 
(^ 1000 A us. 1 — 2A), leading to a shorter exposure time with a trade-off in acctiracy 
of the measured redshifts. 

Considering the new generation of 8m class telescopes, the planned instrumentation 
on the HST and future, high sensitive radio telescopes, a substantial number of deep 
siirveys at different wavelengths will soon become available. Most of these galaxies will 
be too faint for spectroscopic study and hence, the photometric redshift technique is 
the only practical way for estimating their redshifts. In a recent assessment of different 
photometric redshift techniques, using a redshift-limited spectroscopic survey, it was 
shown that photometric redshifts can be predicted with an accuracy of 0.1 (0.3) for 68% 
(99%) of the sources examined (Hogg et al 1998). Therefore, photometric redshifts could 
provide a powerful tool for statistical studies of evolutionary properties of galaxies and 
in particular of faint galaxies for which spectroscopic data are difficult to obtain. 

The most important step in any study concerning photometric redshift measurement, 
is the choice of the template Spectral Energy Distributions (SEDs) for different popu- 
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lations of galaxies, with which the observed SEDs should be compared. There are two 
general ways for adopting the template SEDs: 

a) . Empirical templates: in this case one uses the mean observed SEDs corresponding 
to different types of galaxies. The problem here is that there are not enough information 
about the observed SEDs for different classes of objects at different redshifts (particularly 
at high redshifts). Therefore, incorporating the spectral evolution of galaxies of different 
types on their template SEDs is difficult and uncertain. 

b) . Synthetic templates: uses model SEDs for different spectral types of galaxies, 
shifted in redshift space, assuming evolutionary population synthesis models. The main 
problem here is to constrain the evolutionary model parameters to produce realistic model 
SEDs (for different types) as a function of redshift. In particular, the effect of dust at 
high redshifts (specially in star forming galaxies) is not known. 

To avoid these problems, we introduce a combined approach, producing realistic model 
SEDs based on chemo-photometric Evolutionary Population Synthesis (EPS) models, 
extending from UV to 1 mm in wavelength. The template SEDs here, are consistently 
and simultaneously optimised to; a), produce the observed colours of galaxies at ~ 0; 
b). incorporate chemo-photometric evolution for galaxies of different types, in agreement 
with observations; c). allow treatment of dust contribution and its evolution with redshift, 
consistent with the EPS models; d). include absorption and re-emission of radiation by 
dust and hence, realistic estimates of the far-infrared radiation; e). include correction 
for inter-galactic absorption by Lyman continuum and Lyman forest. The evolutionary 
models and hence, the template SEDs, are constrained by minimising the scatter between 
the photometric and spectroscopic redshifts for a calibrating sample of galaxies with known 
spectroscopic data. 

The main advantage of this technique over the previous works is that it simultane- 
ously and self-consistently allows for the treatment of both the photometric and chemical 
evolution of individual galaxies with time. Moreover, since the synthetic template SEDs 
cover the range from UV to sub- mm wavelengths, one could consistently use the optimised 
SEDs to estimate contributions from individual galaxies to the far-infrared and sub-mm 
wavelengths. Also, the effect of dust and its evolution with redshift is self-consistently 
accounted for in the template SEDs and optimised to produce the observations. This 
is crucial in any photometric redshift technique if it is to be applied on high redshift, 
star-forming galaxies (Meuer et al 1997; Cimatti et al 1998). 

The new photometric redshift technique is outlined in the next section. In section 3, 
the evolutionary population synthesis models are briefly discussed. Section 4 presents the 
calibration sample. This is followed by the optimised template SEDs in section 5. The 
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uncertainties in the photometric redshifts and spectral types are explored in section 6. 
The conclusions are presented in section 7. 

2. Photometric Redshift Technique 

For a given galaxy, the photometric redshift can be estimated by comparing its observed 
SED with a set of template SEDs, corresponding to different morphological types and 
shifted to different redshifts, accounting for galaxy evolution with look-back time. The 
redshift and spectral type associated with the template SED closest to the observed SED 
will then be assigned to that galaxy. The comparison is carried out by minimising the 
function 

n 
i=l 

where the summation, i, is over the passbands with F^f^^ and -F^empiate being, respec- 
tively, the observed and template fluxes at a given passband, n, cr' the uncertainty in 
the observed flux and a the normalisation, estimated in two different ways. First, it is 
estimated by forcing the template SED to have the same total energy (calculated by 
integrating the SED over the observed wavelength range) as the observed SED. Second, 
the normalisation is used as a free parameter in the fit in minimising the x^- The two 
methods give very close results. The comparison between the observed and model SEDs, 
over the entire range of observed SEDs, will constrain the three free parameters (i.e. 
redshift, spectral type and normalisation). However, due to changes with redshift in the 
properties of galaxies, we also need to allow for evolutionary effects on the model SEDs. 
This introduces more free parameters which will be constrained using observed data, as 
discussed in the following sections. 

In order to construct the template SEDs needed for photometric redshift measure- 
ments, we use Evolutionary Population Synthesis (EPS) models, required to estimate 
changes in properties of galaxies with redshift. These models allow us to develope the 
template SEDs as a function of redshift for different types of galaxies. In the next section, 
the EPS models used to derive the template SEDs for different types of galaxies, and 
their underlying assumptions are discussed. The models are constrained, using observed 
local SEDs for different types of galaxies and by optimising the estimated photometric 
redshifts to produce their spectroscopic counterparts for a calibrating sample of HDF 
galaxies with known spectroscopic redshifts. Details of the procedure adopted in this 
study are summarised in the flow chart in Figure 1. 
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3. The Evolutionary Population Synthesis Models 

Synthetic SEDs are constructed following a self consistent procedure, incorporating stel- 
lar emission, internal extinction and re-emission by dust. These chemo-photometric EPS 
models, providing the SEDs extending over four decades in wavelength (ie. from O.OS/im 
to 1 mm)- (Mazzei et al 1992, hereafter MXD92), allow us to investigate the local proper- 
ties and the evolution with galactic age of the SEDs of different galaxy types. In particu- 
lar, we use realistic EPS models for disc galaxies (MXD92), early-type galaxies (Mazzei, 
De Zotti and Xu (1994), hereafter MDX94; Mazzei and De Zotti (1996)) and the far- 
infrared luminous starburst population (Mazzei, Curir and Bonoli (1995)). 

In order to match the photometric properties of the galaxy types, expected to dom- 
inate deep surveys, we consider 4 different templates, consisting of; elliptical, spiral, 
irregular and starburst. Details of the SED models for these types and their evolution 
with cosmic time are given in the above references. In this section, an overview of the 
main parameters (to be optimised using observed SEDs) will be presented. 

The ESP models, used to generate synthetic SEDs for different types of galaxies, have 
been linked to their chemical evolution so that the increased metallicity of successive stel- 
lar generations is taken into account. We adopt a parametrization for the time evolution 
of the star-formation rate as 

m = Mg{t) M^/yr, (1) 

where fg is the fractional mass of gas which takes part in the star formation {fg = 
^gas I ^^gai) cind t/'q is thc initial star-formation rate ( i.e. SFR, at z — z form 

; where 

Zform is the formation redshift). Wo assume that Initially fg = 1 with m.gai = 10^^ Mq. 

Thc synthetic spectra for thc stellar generations with different mctalllcltles which 
contribute to thc galaxy SED are estimated at any given time, using the re- 
cent set of Isochroncs constructed by BcrtcUl ct al. (1994). These Incorporate the 
results of new stellar evolutionary calculations, based on six metallicity values: 
0.0004,0.001,0.004,0.008,0.02 and 0.05. The Isochroncs have a fine coverage of masses 
and ages and Include almost all the evolutionary phases from the main sequence to the 
stage of planetary nebulae ejection or carbon Ignition, as appropriate for a given Ini- 
tial mass. This allows us to link the photometric and chemical evolution In galaxies. 
Further, at any given time and wavelength, the models Incorporate both the average 
correction due to Internal extinction by the enriched ISM, assuming a dust to gas ratio 
proportional to the gas metallicity and the dust re-radlatlon from 3 /im up to 1 mm. We 
assume that gas and stars have the same distribution. In particular, for spiral, starburst 
and Irregular galaxies we follow the same approach as described In MXD92 (I.e. an ex- 
ponential function of the galacto-centrlc radius- see their eq.(ll)) whereas for ellipticals 
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we follow MDX94 (i.e a density profile given by the King (1966) model). Results are 
not strongly dependent on the adopted geometry. The evolution of the optical depth, 
T, follows directly from that of the gas metallicity, Zg^s and fg since T(t) oc Zgas{tyfg 
( as a consequence of the assumption of a gas-to-dust ratio being proportional to the 
gas metallicity; MXD92, MDX94). This results in very different histories for the effective 
optical depths (where Lpm/Lboi = 1 — exp{—Teff)) in the early phases of the evolution 
of ellipticals in the model, which is based on the adopted IMF parameters. 

The IMF and the SFR are the main parameters, controling the e-folding star forma- 
tion time scale, the chemical evolution and, as a direct consequence, the optical depth of 
the system. Thus, the far-infrared data, where available, are expected to constrain the 
galaxy metallicity. Therefore, by comparing the local observed SEDs of different galaxy 
types with models, we constrain both Vo and the IMF. The observed SEDs for nearby 
spiral and elliptical galaxies with known metallicities and far-IR (IRAS) observations 
are used to constrain the IMF parameters so that these models successfully reproduce 
the local chemo-photometric properties of both these types over four decades in wave- 
length. However, the parameter which is most sensitive to the local properties (i.e. the 
gas content and optical colours) of galaxies along the Hubble sequence, is the initial 
star formation rate, tpo (Sandage 1986), which will be constrained for different types of 
galaxies, using the observed SEDs (section 5.1). 

4. The Calibration Sample 

To further constrain the evolutionary parameters in these models, in particular n and 
Zform, a calibration sample of HDF galaxies with available spectroscopic rcdshifts is 
compiled. The calibration sample contains galaxies bright enough to allow spectroscopic 
redshift measurements, consisting of both objects with UBVI detections and the UV 
'drop outs' (i.e. high redshift objects). The advantage of using the HDF galaxies as the 
calibration sample here is two- fold: (a). HDF galaxies cover the spectral range required 
for photometric redshift measurement, containing galaxies in the range < z < 3.5; (b). 
they are selected using a uniform criteria, for both nearby and distant galaxies, making 
this an unbiased calibration sample. 

The calibrating galaxies are individually inspected and for galaxies with close neigh- 
bours, a smaller aperture is adopted to avoid contamination of their light by nearby 
objects. For the objects with no UV detection (i.e. UV 'drop outs'), a U-band magnitude 
of 28.01 mag. was assumed for the photometric redshift measurement. Although this is 
likely to introduce a bias due to a colour-magnitude relation (i.e. fainter galaxies are 
bluer), we do not expect it to be significant. This is examined by exploring the range 
28.0 < U < 30.0 for each of the galaxies in the UV- drop out sample and estimating 
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their respective photometric redshift. Using models with the same parameters, we find, 
on average, only a small sensitivity of the estimated photometric redshifts (< 3a) to 
changes in the UV magnitudes in the above range. 

The calibrating sample, consisting of 73 galaxies, is listed in Table 1 together with 
their UBVI photometry and spectroscopic redshifts. The reliability of the spectroscopic 
redshifts and the photometric accuracy of individual galaxies are discussed in the foot- 
note to this table. The calibration sample in Table 1 was selected to be in the HDF area, 
to allow accurate UBVI photometry and to have unambigious spectroscopic redshifts. 
Moreover, objects with non-stellar sources of energy (i.e. gravitationally lensed candi- 
dates; Zepf et al (1996)) are not included. The magnitudes are in the AB system and are 
measured over an aperture of 3 arcsec diameter (unless stated otherwise in the footnotes 
to Table 1). 

The EPS models for different types of galaxies are therefore further constrained by 
minimizing the rms scatter between the photometric redshifts, predicted by our models, 
and their spectroscopic counterpart, using the calibration sample in Table 1. Thus, the 
final EPS models for different types of galaxies, which we define as templates (see below), 
have local {z = 0) SEDs which match the local, observed, SEDs of the corresponding 
galaxy types and evolutionary properties constrained, using the calibration sample here. 

5. The Template SEDs 
5.1. EPS Model Parameters 

A large number of EPS models with different input parameters (IMF, tpQ, n and Zforrn) 
are developed, each aceounting for both the local properties and evolutionary behaviour 
of the 4 types of galaxies considered here. For a given spectral type, the model parameters 
are normalised at z = by fitting them to the local observed SED of their respective type. 
The evolutionary behaviour of the EPS models are then constrained by estimating the 
photometric redshifts to our calibrating sample in Table 1, considering template SEDs 
with different parameters (i.e. IMF and z/orm) and allowing for Lyman continuum and 
Lyman forest absorption for z > 2 galaxies. The template SEDs corresponding to the 
evolutionary parameters which give the closest agreement between the photometric and 
spectroscopic redshifts are then adopted (see below). The sensitivity of the final results 
(i.e. photometric redshifts) to the input parameters in the EPS models is studied in the 
next section while, details of the final templates for individual types, which best satisfy 
the above requirements, are summarised below: 

a). Ellipticals. The synthetic SEDs, representing elliptical galaxies, are constructed 
taking n ~ 0.5 and ipo = 100 MQ/yr. This gives an e-folding star formation time scale of 1 
Gyr and reproduces the observed local SEDs for the ellipticals (Figure 2.1a). A formation 
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redshift of 5 is estimated, corresponding to an age of 13 Gyrs. These models predict a 
noticeable extinction by dust (As < 4.6 mag) in the first stages of their evolution which 
are characterized by intense star formation activity, making them powerful far-infrared 
sources (Fig 2.2a) (see MDX94 for more details). 

b) . Spirals. The models used to produce synthetic SEDs for spiral galaxies have n = 2 
and tjjo = IOMq/yt. They are consistent with an e- folding star formation time scale of 
10 Gyrs, corresponding to a formation redshift of 2 (i.e. an age of 10-11 Gyrs). This gives 
a local value of Lfir/Lbol = 0.3, very similar to that of our own Galaxy. This model 
re-produces the local SED of NGC 3627 out to the far-infrared wavelengths (Figure 2.1c). 
Spiral galaxies are slowly evolving with time {Lfir/Lbol < 0.1 beyond z = 1), resulting 
a smooth evolution for their SED (Figure 2.2c)- (see MXD92 for more details). 

c) . Starbursts. The template SEDs for the starburst population are produced taking 
n = —1, tpo = 2MQ/yr and a formation redshift of 5. This template does not produce 
powerful far-infrared emission at any redshift and hence, its local SED (Figure 2.1b) is 
different from that of local luminous far-infrared starbursts (ie. M82 and Arp 220). The 
SFR in this model is a gradual process with a smooth time scale, leading to formation of 
very blue, metal poor systems at 2; ~ (Figure 2.1b) and a blue, dust-free {Lfir/Lbol < 
0.05) system at > 1 (Figure 2.2b). This mimics a scenario involving frequent but short 
bursts of star formation, which use a small fraction of gas in these systems (strong bursts 
of star formation rapidly exhaust the gas, leading to ellitpical like systems). The starburst 
templates here represent the average evolutionary behaviour expected for this population 
of galaxies. 

d) . Irregulars. The Irregular template has been produced with the same recipe as the 
starburst but with a formation redshift of Zform = 1 (ie. an age of 0.8-0.9 Gyr). The 
local SED of the irregular galax;y NGC 4449 (Kennicutt, 1992) is compared in Figure 
2. Id with our local templates for both irregular and starburst galaxies. These templates 
are also compared at 2; = 0.975 (Figure 2. 2d), showing a significant difference beyond 
A lyum. 

The model SEDs at z = for the four galaxy types discussed above, agree well with 
the local observed SEDs. as shown in Figures 2.1a-2.1d. The effects of Lyrnari break and 
Lyman forest opacities are included to the template SEDs, using the relations A{M) vs. 
z for different wavelengths, given in Madau et al (1996). These relations were fitted to 
parametric forms, which were then used to estimate the respective correction (due to 
absorption by inter galactic medium) to the SEDs at any given redshift. The correction 
due to Lyman break and Lyman forest absorption ranges from A{Muv) = 0.25 mag. at 
z = 2 to A{Muv) ~ 1 mag. at 2; = 2.5. 
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5.2. Sensitivity of the Template SEDs on the Model Parameters 

In this section we study the dependence of the results (i.e. the photometric redshifts) to 
the model parameters which most strongly affect the final template SEDs and hence, the 
predicted photometric redshifts. These parameters consist of the shape of the IMF and 
its lower mass limit, the total number of templates (i.e. spectral types) and the formation 
redshift for each galaxy type. New template SEDs are generated corresponding to EPS 
models for elliptical, spiral, starburst and irregular galaxies, using the parameters listed 
in Table 2, with the rest of the parameters taken to be the same as discussed in section 5.1. 
For each set of the new templates, the photometric redshifts are estimated for galaxies in 
the calibrating sample (Table 1) with the rms scatter in the quantity {zphot — -Zspec)/(1 + 
Zspec) (ie. between the photometric and spectroscopic redshifts) calculated and listed in 
Table 2. The templates corresponding to the model which gives the smallest rms estimate 
in Table 2 (i.e. model 4) is then adopted. The photometric redshifts estimated for the 
calibration sample, using the adopted template SEDs (model 4 in Table 2), are listed in 
Table 1 and compared with their spectrosopic counterparts in Figure 3. The rms scatter 
of 0.11 here is taken as the uncertainty in the photometric redshift estimates in the 
range < ^; < 3.5. The uncertainties in correcting for intergalactic absorption at 2; > 2 
indirectly affect the optimisation of the EPS model parameters in this section. To explore 
this, we constrained the calibration sample only to galaxies with z <2 (which are much 
less affected by the IGM absorption) and estimated the EPS model parameters so that to 
minimise the rms scatter in Figure 3. No change is found in the EPS model paremeters. 
Also, we explored the sensitivity of the rms scatter in Figure 3 to the number of template 
SEDs used (ie. including templates for different sub-classes of spirals), taking the number 
of SEDs as a free parameter. This did not reduce the optimum rms scatter, derived using 
the four templates. 

Considering the results in Table 2, it seems that one set of models (Salpeter IMF with 
TO; = 0.01), give a considerably better fit to observations. This was extensively tested 
by exploring the parameter space, consisting of the IMF shape and its mass limits, total 
number of templates and formation redshifts. It was found that this is not due to sampling 
a particular region of the parameter space or the number of templates, and is indeed, 
a real effect. A similar study, using a different set of EPS models and the optimised 
parameters in Table 2, would be extremely valuable. 

The galaxies in the calibration sample in Table 1 also have near-IR data (Fernandez- 
Soto et al. 1998). The above procedure was repeated, using the combined optical and 
near-IR magnitudes (UBVIJK) for the calibrating sample. This did not change the opti- 
mized EPS model parameters in section 5.1 and Table 2. As an independent test of the 
template SEDs here, we include the near-IR magnitudes to the observed SEDs of the cal- 
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ibrating sample in Table 1 and estimate their photometric redshifts, using the template 
SEDs predicted in section 5.1 (model 4 in Table 2). Compared to their spectroscopic 
counterparts, an rms scatter of 0.13 is found, in agreement with 0.11 from Figure 3. 

The amount of dust and its evolution with redshift is an important characteristic of the 
EPS models in this study. This, at any time, is computed self-consistently, accounting 
for both the SFR and the IMF parameters. By extending the IMF to low m; values 
{mi = 0.01 Mq), the gas depletion rate becomes faster, reducing the dust enrichment 
rate. This leads to higher optical depth in the early evolutionary phase of our elliptical 
models (templates)- (see also Mazzei & DeZotti 1996 for more details). The UV extinction 
{Au) corresponding to different IMF and mi values is estimated for both elliptical and 
spiral templates at different redshifts and are listed in Table 3. By changing the shape 
of the IMF and its lower mass limit, the UV extinction for ellipticals at ^ ^ 2 changes 
in the range ~ 1 — 4 mag. Considering the optimised model in Table 2 (model 4), the 
extinction in spirals is significantly larger than in ellipticals at 2: ~ while, at higher 
redshifts {z ~ 1.5), ellipticals are obscured more than the spirals. Using our optimised 
model in section 5.1, the template SEDs at 2 ~ 1.5 are compared in Figure 2.2 with their 
counterparts at 2; ^ 0. These show a significant dust contribution to the elliptical SEDs 
at 2; ~ 1.5, as indicated from the peak at the far-IR wavelengths. 

6. Uncertainties in the Photometric Redshifts and Spectral Types 

A potential source of uncertainty in estimating the photometric redshifts and spectral 
types of galaxies is the possibility that the models might be degenerate (i.e. different 
synthetic SEDs, corresponding to different redshifts and galaxy types producing the 
same result). Furthermore, the photometric errors in the observed SEDs are likely to 
affect the final estimate of both the photometric redshifts and the spectral types of their 
respective galaxy. Also, due to the relative similarity of the model SEDs for star burst 
and irregular galaxies (Figure 2), the accuracy with which the spectral types for these 
systems are predicted, needs to be established. 

To investigate the above problems, a Monte Carlo simulation is performed. A simu- 
lated catalogue is generated to resemble the observed HDF survey, with UBVI magni- 
tudes, known redshifts (zinput) and spectral types. The galaxies are randomly selected 
to have SEDs similar to the synthetic SEDs for the four types of galaxies (elliptical, spi- 
ral, starburst and irregular) considered in section 5.1, shifted in redshift space. Random 
Gaussian noise, resembling photometric errors, are then added to the simulated SEDs. 
The simulated catalogue has a magnitude limit of 1=28 mag and an apparent magnitude 
distribution similar to the observed HDF survey. 
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The photometric redshift code is used to predict the redshifts (zoutput) and spectral 
types of individual galaxies in the simulated catalogue and to compare them with their 
input values. Two sets of simulations is carried out, assuming different observational 
errors (i.e. a in the Gaussian noise distribution). The difference between the input (sim- 
ulated) and the output (predicted) redshifts is presented in the histogram in Figure 4, 
with the spectral types compared in Table 4. 

The log{zoutput/zinput) distribution (Figure 4), including all the four types of galaxies, 
shows a distinct peak at zero, indicating that the redshifts for the simulated galaxies are 
well re-produced within Az 0.11. The objects, located at the tails of the distribution 
in Figure 4, are all galaxies classified as spirals, irregulars or starbursts for which, the 
degeneracy of their SEDs appears to be more serious (see below). When increasing the 
photometric errors in the simulated SEDs, the agreement between the simulated (input) 
and predicted (output) redshifts decreases but the distribution still peaks around zero 
(dotted line in Figure 4). 

The spectral types of the simulated galaxies are compared with the predicted types 
in Table 4. This gives the ratio of galaxies of a given spectral type, which are correctly 
classified in the simulation (i.e. galaxies for which their spectral type in the input cat- 
alogue were successfully re-produced), to the total number of galaxies of the same type 
in the input catalogue. In both simulations, we re-produce the spectral types for all the 
ellipticals in the input catalogue with no mis-identifications. However, due to the rela- 
tive similarity of the UV-to-optical part of the synthetic SEDs for the spirals, irregulars 
and starburst galaxies (Figure 2), we can recover respectively 79%, 85% and 71% of the 
spectral types of galaxies in the input catalogue. 

The conclusion from the simulation here is that the photometric redshifts are, on 
average, well produced (within the expected accuracy) and are not affected by the de- 
generacy in the template SEDs. Furthermore, the spectroscopic type classification for 
ellipticals is reliable with 100% re-produced (i.e. no mis-identifications) from the input 
catalogue while, for the spirals, irregulars and starbursts, there is a slight degeneracy in 
predicting their spectral types. 

7. Summary and Conclusion 

A now technique is developed for estimating the photometric redshifts to galaxies. The 
advantage of this over the previous methods is that it allows a self-cosistcnt treatment 
of both photometric and chemical evolution of the template SEDs for different types of 
galaxies, in agreement with observations. For this reason, the present technique is partic- 
ularly useful for application to high-redshift galaxies. The template SEDs are constrained 
to simultaneously satisfy the following criteria: a), produce the observed colours of galax- 
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ies at z = 0; b). incorporate the spectral evolution for galaxies of different types; c). allow 
a self-consistent treatment of the dust contribution and its evolution with redshift; d). 
include correction for inter-galactic absorption by Lyman continuum and Lyman forest. 
The model parameters are then constrained, using a a calibrating sample of HDF galax- 
ies. Using only four optical passbands (UVRI), the photometric redshifts to galaxies 
can be estimated with an rms accuracy of 0.11. Including near-IR data (UVRIJK), gives 
a similar rms scatter in the photometric redshifts from this method. The sensitivity of 
the results to different EPS model parameters are explored and the degeneracy of the 
photometric redshifts and spectral types for galaxies are examined using a simulated 
catalogue. 
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Fig. 1. Flow chart summarizing the photometric redshift technique. 

Fig. 2. Fig. 2.1). The synthetic Spectral Energy Distributions (SEDs) are predicted and 
constrained to fit the observed data at z ~ 0, as explained in the text (left panel). Details 
of each panel and the source of the observational data for each type is given below. 

a) . Elliptical galaxies: open squares (Burstein et al. (1988)); filled triangles (Schild and 
Oke (1971)); filled circles (Oke and Sandage (1968)); asterix (Kennicutt (1992)); the filled 
circles at FIR wavelengths correspond to the average local FIR SED for these galaxies 
(Mazzei and DeZotti (1994b)). 

b) . Starburst galaxies: this corresponds to the observed SEDs of NGC5996 (MK691). 
The observed data arc taken from: filled circles (Kennicutt (1992)); asterix in the shorter 
wavelength region (lUE data from Kinney et al. (1993)); asterix in the near-IR region 
(Balzano (1983)); FIR data (IRAS Catalugue Version 2 (Fullmer, L. and Lonsdale, C. 
1989)); the lUE data are measured over 20" x 10" aperture and are shifted vertically 
by a factor of 1.6 to normalise to the optical SED; near-IR data, measured over 10.3" 
aperture has been shifted by a factor of 4.5. 

c) . Spiral galaxies: this corresponds to the NGC3627 galaxy. The observed data are from 
Kennicutt (1992) and Rice et al. (1988). 

d) . The local observed SED of IRR galaxy NGC 4449 (Kennicutt, 1992), filled circles, is 
compared with the synthethic SED of our Irr template at z = (continuous line) and 
that for the starburst galaxies (dashed line). 

Fig. 2.2a-c. The local synthetic SEDs for different types of galaxies from Figures 2.1a-c 
(solid lines) are compared with their counterparts at 2; ~ 1.75 (dashed lines). 
The SEDs for irregulars (solid line) and starbursts (dashed line), both at 2; = 0.975, are 
compared in Figure 2. 2d. 

Fig. 3. Comparison between the spectroscopic redshifts with the photometric values 
estimated here for 73 HDF galaxies from Table 1. The spectroscopic redshifts are taken 
from Cohen et al (1996); Steidel et al (1996); Lowenthal et al (1997); Dickinson (1998). 
The UV drop-out objects for which a U-band limiting magnitude of 28.01 is used to 
estimate the photometric redshifts, are also included. The crosses are HDF3646-I-1408 
and HDF3659-I-1222 galaxies which have uncertain published spectroscopic redshifts of 
2; = 0.13 and 0.47 respectively. The lines have a slope of unity. The photometric redshifts 
are based on model 4 in Table 2. The rms scatter is estimated in {zphot — z spec)/ + z spec) 
and corresponds to 0.11 with the dashed lines corresponding to ±3a error. 
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Fig. 4. Histogram of log{zoutput/ Zmput) for the simulated UV selected HDF catalogue. 
The continous line is the result for the photometric errors (Af/, AB, AF, A7) corre- 
sponding to (0.15,0.10,0.05,0.05) and the dotted line corresponds to (0.25,0.15,0.10,0.10). 
Both the histograms peak around zero, implying that the redshifts in the simulated cat- 
alogue can be produced to within A.z = 0.11. 
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Table 1. The Calibrating Sample 





RA(J2000) 


Dec(J2000) 


I 


R 


B 


U 


^spec 


^phot 


HDF3646+1143 


12 


36 


46.12 


62 


11 


43.38 


21.47 


22.44 


22.92 


23.25 


1.012 


0.775 


HDF3643+1143 


12 


36 


43.72 


62 


11 


43.95 


20.88 


22.49 


24.62 


24.62 


0.764 


0.700 


HDF3646+1146 


12 


36 


46.84 


62 


11 


46.09 


23.50 


24.17 


24.24 


23.89 


1.059 


0.950 


HDF3648+1147 


12 


36 


48.31 


62 


11 


47.47 


25.12 


25.28 


26.04 


28.01 


2.980 


3.000 


HDF3646+1152 


12 


36 


46.44 


62 


11 


52.31 


21.91 


22.99 


24.70 


24.58 


0.504 


0.425 


HDF3643+1152 


12 


36 


43.32 


62 


11 


52.63 


22.91 


23.33 


23.42 


23.61 


1.242 


1.600 


HDF3649+1156 


12 


36 


49.26 


62 


11 


56.08 


23.48 


23.92 


23.80 


23.48 


0.961 


0.975 


HDF3641+1201 


12 


36 


41.51 


62 


12 


1.41 


24.64 


25.03 


25.83 


25.47 


0.483 


0.650 


HDF3645+1202 


12 


36 


45.86 


62 


12 


2.63 


23.66 


24.44 


24.83 


24.26 


0.679 


0.900 


HDF3640+1203 


12 


36 


40.70 


62 


12 


3.84 


23.10 


23.73 


23.97 


24.44 


1.010 


1.350 


HDF3641+1206 


12 


36 


41.85 


62 


12 


6.59 


21.01 


21.64 


22.57 


23.19 


0.432 


0.600 


HDF3651+1210 


12 


36 


51.94 


62 


12 


10.69 


22.92 


23.44 


24.17 


24.81 


0.456 


0.525 


HDF3648+1215 


12 


36 


48.17 


62 


12 


15.04 


22.56 


23.40 


23.75 


23.58 


0.960 


0.875 


HDF3642+1217 


12 


36 


42.78 


62 


12 


17.48 


20.73 


21.30 


22.20 


22.88 


0.454 


0.550 


HDF3643+1219 


12 


36 


43.53 


62 


12 


19.33 


22.59 


23.42 


23.92 


24.71 


0.750 


0.725 


HDF3652+1220 


12 


36 


52.60 


62 


12 


20.77 


23.15 


23.63 


24.37 


25.12 


0.401 


0.480 


HDF3649+1221 


12 


36 


49.44 


62 


12 


21.17 


24.04 


24.50 


24.65 


24.45 


0.961 


0.950 


HDF3651+1221 


12 


36 


51.62 


62 


12 


21.30 


21.54 


22.23 


23.32 


24.61 


0.299 


0.460 


HDF3659+1222'* 


12 


36 


59.45 


62 


12 


22.20 


23.67 


24.30 


24.68 


24.18 


0.470 


0.925 


HDF3648+1222 


12 


36 


48.95 


62 


12 


22.38 


22.43 


23.21 


23.53 


23.50 


0.953 


0.875 


HDF3658+1222 


12 


36 


58.56 


62 


12 


22.76 


23.34 


24.03 


24.40 


23.90 


0.681 


0.925 


HDF3657+1226 


12 


36 


57.11 


62 


12 


26.68 


22.45 


23.03 


23.64 


23.42 


0.561 


0.700 


HDF3651+1228 


12 


36 


51.26 


62 


12 


28.47 


26.33 


25.94 


26.41 


27.30 


2.775 


2.800 


HDF3644+1228 


12 


36 


44.55 


62 


12 


28.51 


23.59 


23.96 


24.10 


25.37 


2.268 


2.050 


HDF3647+1231 


12 


36 


47.30 


62 


12 


31.76 


22.82 


23.23 


23.99 


23.95 


0.421 


0.550 


HDF3653+1235 


12 


36 


53.37 


62 


12 


35.41 


22.81 


23.48 


24.21 


26.13 


0.559 


0.430 


HDF3700+1235 


12 


37 


0.48 


62 


12 


35.77 


21.40 


22.54 


24.11 


24.30 


0.562 


0.440 


HDF3646+1237 


12 


36 


46.96 


62 


12 


37.95 


21.08 


21.50 


22.26 


22.91 


0.321 


0.175 


HDF3641+1239 


12 


36 


41.63 


62 


12 


39.82 


25.48 


25.51 


25.53 


28.13 


2.591 


2.300 


HDF3650+1240 


12 


36 


50.12 


62 


12 


40.94 


20.76 


21.31 


22.14 


22.76 


0.474 


0.550 


HDF3644+1241 




OD 


/I/I ns 

44. Uo 


DZ 


1 




23.16 


23.91 


24.16 


24.21 


0.873 


0.900 


HDF3655+1246 


12 


36 


55.49 


62 


12 


46.35 


21.95 


22.93 


23.58 


23.79 


0.790 


0.750 


HDF3650+1246 


12 


36 


50.18 


62 


12 


46.92 


21.45 


22.89 


24.58 


26.60 


0.678 


0.625 


HDF3644+1248 


12 


36 


44.11 


62 


12 


48.97 


21.49 


21.98 


22.55 


22.94 


0.558 


0.625 


HDF3649+1249* 


12 


36 


49.74 


62 


12 


49.94 


25.17 


25.44 


26.62 


28.01 


3.233 


3.200 


HDF3655+1250 


12 


36 


55.57 


62 


12 


50.01 


23.12 


23.84 


24.06 


23.92 


0.950 


0.925 


HDF3658+1253 


12 


36 


58.66 


62 


12 


53.26 


21.50 


21.95 


22.76 


23.17 


0.320 


0.470 


HDF3656+1253 


12 


36 


56.66 


62 


12 


53.56 


23.18 


23.99 


24.28 


24.03 


1.231 


0.900 


HDF3647+1253 


12 


36 


47.47 


62 


12 


53.68 


23.74 


24.36 


25.12 


25.03 


0.681 


0.700 
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Table 2. Sensitivity of the rms scatter between the photometric and spectroscopic red- 
shifts on the EPS model parameters (i.e. template SEDs) 



Model 


IMF 




Number of 






orm 




rms 








Templates 


E 


Sp 


Sb 


Irr 




1 


Salpeter 


0.01 


4 


5 


5 


5 


1 


0.18 


2 




0.01 


3 


5 


5 


5 




0.26 


3 




0.01 


3 


5 


2 


5 


_ 


0.24 


4 




0.01 


4 


5 


2 


5 


1 


0.11 


5 


Salpeter 


0.10 


4 


5 


5 


5 


1 


0.27 


6 




0.10 


3 


5 


5 


5 




0.38 


7 


5) 


0.10 


3 


5 


2 


5 




0.37 


8 




0.10 


4 


5 


2 


5 


1 


0.22 


9 


Scalo 


0.10 


4 


5 


5 


5 


1 


0.28 


10 




0.10 


3 


5 


5 


5 




0.52 


11 


5) 


0.10 


3 


5 


2 


5 




0.45 


12 


75 


0.10 


4 


5 


2 


5 


1 


0.25 



Table 3. Changes in UV extinction (Au) with redshift, corresponding to EPS models 
for Elliptical and Spiral galaxies. 

Salpeter IMF Salpeter IMF Scalo IMF 





mi = 


0.01 M0 


mi = 


0.1 Mq 


m.i = 


0.1 M( 


redshift 


AuiE) 


AuiSp) 


AuiE) 


Au (Sp) 


AuiE) 


Aui: 





.02 


0.78 


.02 


2.18 


.02 


2.45 


0.5 


.44 


0.35 


.34 


1.22 


.29 


1.44 


1 


0.67 


0.15 


.56 


0.57 


.43 


0.70 


1.5 


2.32 


0.04 


1.48 


0.16 


.83 


0.20 


2. 


4.23 




2.23 




1.17 




2.5 


4.79 




2.28 




1.25 
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Table 4. Results from the simulation of spectral types. Fraction of the galaxies of 
different types, correctly classified between the input and output simulated catalogues 
are listed. 

Photometric Errors Percentage^ 
A{U) A{B) A{V) A{I) E Sp Sb Irr 

simulation 1 0.15 0.10 0.05 0.05 100% 79% 71% 85% 

simulation 2 0.25 0.15 0.10 0.10 100% 66% 49% 62% 

^ Ratio of the number of galaxies of a given spectral type in the output catalogue which are 
correctly classified, to the total number of galaxies of the same spectral type in the input 
catalogue. 



Model SEDs (E, Sp, Sb, Irr) 

[IMF, t/jojg^n, Zgas, Zform] 




Optimized to fit local 
observed SEDs [/^(-^o)? 







Optimize model parameters 
by minimizing rms 
between ZpUot and Zgp^c 



fg — fTT-gas / ^gal 



Estimate Zphot for 
galaxies with known Zgpec 
for given model parameters 
(IMF, 

Zform etc) 



^phot ^spec 
^~\~^spec 



Final Template SEDs 
for E, Sp, Sb, Irr 
rms=0.11 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-J \* ^ 

- °° r f"^ - 

_ a 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

/ " - 
/ \ - 
X / ^ - ^ \ 

r \ 
z=0 

z=1 .75 _ _ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
_ b ^^^""X 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

\ \ 
\ \ 

\ \ ~ 

\ \ 
\ \ 

\ \ - 

z=1.75 _ _ \ 

\ \ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

_ c Y 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

\ ~ 
\ 

"X / V\ 
\y \\ 
\\ 

\\ ~ 
\\ 

\\ _ 
z=0 ^T^^^^ _ 
z=1.75 _ _ - ^ Y 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 L 1 1 



1 -0.5 0.5 1 1.5 2 -1 -0.5 0.5 1 1.5 2 

log (A yu.m) log (A /i.m) 



2 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




-1 -0.5 0.5 1 1.5 2 

log (A /xm) 



150 



125 



100 



75 



50 



25 





-0.25 -0.2 -0.15 -0.1 



-0.05 0.05 0.1 

log(z„,tput/Zinpul) 



0.15 0.2 0.25 



